***** INEP_LOGIT AILP
*importar dataset #####alterar diretorio#####
import excel "C:\Users\User\Navus\datasetAILP.xlsx", sheet("Planilha1") firstrow

*caracteristica do dataset
desc

*Distribuicao de frequencia das variaveis qualitativas
tab grupolp

*Regressao logistica binominal
logit grupolp nivel_socio_economico mird ied_ai icg had_f14 afd_ai atu_f14, level(95)
	
*Salvar estimativas dos parametros do modelo completo
lrtest, saving(0)

*Regressao logistica stepwise forward
stepwise, pr(0.10) pe(0.05) forward: logit grupolp nivel_socio_economico mird ied_ai icg had_f14 afd_ai atu_f14, level(95)

*Odds ratio
logit, or

*likelihood-ratio test (teste de razão de verossimilhanca)
lrtest

*AIC e BIC
estat ic

*Teste de Hosmer-Lemeshow
estat gof, group(10) table

*Sensibilidade do modelo
estat class, cutoff(0.50)

*Curva de sensibilidade
lsens

*Curva ROC
lroc

*DFBETA
predict dbeta, dbeta
scatter dbeta id_escola, mlab(id_escola) 
sum dbeta

*VIF
reg grupolp nivel_socio_economico mird icg, level(95)
vif

*Reestimacao de Jackknife
jackknife _b, level(95) : logit grupolp nivel_socio_economico mird icg, level(95)


***** INEP_LOGIT AIMT
*importar dataset #####alterar diretorio#####
import excel "C:\Users\User\Navus\datasetAIMT.xlsx", sheet("Planilha1") firstrow


*caracteristica do dataset
desc

*Distribuicao de frequencia das variaveis qualitativas
tab grupomt

*Regressao logistica binominal
logit grupomt nivel_socio_economico mird ied_ai icg had_f14 afd_ai atu_f14, level(95)

*Salvar estimativas dos parametros do modelo completo
lrtest, saving(0)

*Regressao logistica stepwise forward
stepwise, pr(0.10) pe(0.05) forward: logit grupomt nivel_socio_economico mird ied_ai icg had_f14 afd_ai atu_f14, level(95)
*stepwise, pr(0.10) pe(0.05) forward: logit grupomt nivel_socio_economico mird ied_ai icg had_f14 afd_ai atu_f14, level(95) noconstant

*Odds ratio
logit, or

*likelihood-ratio test (teste de razão de verossimilhanca)
lrtest

*AIC e BIC
estat ic

*Teste de Hosmer-Lemeshow
estat gof, group(10) table

*Sensibilidade do modelo
estat class, cutoff(0.50)

*Curva de sensibilidade
lsens

*Curva ROC
lroc

*DFBETA
predict dbeta, dbeta
scatter dbeta id_escola, mlab(id_escola) 
sum dbeta

*VIF
reg grupomt nivel_socio_economico mird icg atu_f14, level(95)
vif

*Reestimacao de Jackknife
jackknife _b, level(95) : logit grupomt nivel_socio_economico mird icg atu_f14, level(95)
